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I.  Summary  of  Research  Highlights 

Our  work  studies  novel  phases  of  matter  potentially  realizable  with  ultracold  atomic  gases 
confined  in  optical  lattices  and  models  the  impact  of  controlled  and  uncontrolled  disorder 
on  these  optical  lattice  phases.  Atoms  (either  bosons  or  fermions)  placed  in  optical  lattices 
can  interact  to  form  strongly  correlated  states.  Our  work  seeks  to  foster  experimental  work 
with  these  states,  which  include,  superfluids,  superconductors,  supersolids,  glasses, 
topological  states,  and  other  quantum  many-body  states  because  they  are  of  fundamental 
importance  in  understanding  the  collective  behavior  of  quantum  particles.  These  states  are 
potentially  impacted  by  spatial  and  phase  fluctuations,  trapping,  laser  noise,  and  other 
sources  of  uncontrolled  disorder.  Controlled  disorder  can  also  be  used  to  drive  the  system 
into  interesting  new  regimes  as  well. 

This  award  has  led  to  significant  accomplishments  over  the  past  three  years  that  pertain 
directly  to  the  goals  of  the  grant: 

•  22  publications  and  manuscripts  (5  in  Physical  Review  Letters) 

•  13  students  and  postdocs  mentored 

•  27  invited  presentations  given 

The  award  has  also  allowed  us  to  establish  direct  connections  with  experimental  groups  (via 
travel  and  shared  projects).  These  groups  include  Peter  Engel’s  group  at  Washington  State 
University,  Ian  Spielman’s  group  at  NIST /JQI/UMD,  and  Brian  DeMarco’s  group  at  the 
University  of  Illinois  at  Urbana-Champaign. 

Our  work  established  and  analyzed  models  that  guide  the  ultracold  atom  community  in 
realizing  novel  states  of  matter  with  optical  lattices.  Our  work  discovered  new  classes  of 
many-body  models  that  meet  goals  of  the  proposal  1)  Our  models  of  controlled  disorder 
show  that  it  can  be  an  asset  to  revealing  many-body  states  2)  Strong  interactions  can  often 
stabilize  intriguing  states  of  matter  against  fluctuations  (particularly,  thermal  fluctuations). 
3)  Topological  quantum  phases  are  robust  against  quantum  and  thermal  fluctuations.  The 
following  briefly  summarizes  some  of  our  major  research  accomplishments.  Section  V 
expands  on  this  summary  as  well. 

•  The  Bose  glass  is  an  intriguing  disordered  state  of  bosons  defined  by  local  order.  To- 
date  it  has  not  been  observed  experimentally.  We  proposed  a  local  measure,  boson 
core  compressibility,  which  uses  double  occupancy  for  identification  of  a  Bose-Glass 
and  other  states  (Publication  18). 

•  A  supersolid,  a  simultaneous  superfluid  and  solid,  might  be  made  possible  via  dipolar 
interactions  between  bosons  in  optical  lattices.  It  has  not  been  observed  to  date. 
Supersolids  are  believed  to  be  rather  delicate,  but  we  discovered  that  controlled 
disorder  enhances  the  strength  of  supersolids  (Publication  19). 

•  Topological  states  of  matter  can  reveal  intriguing  excitations,  Majorana  fermions, 
that  interact  non-locally  to  demonstrate  anyon  braid  statistics.  Unfortunately,  their 
non-local  character  is  not  observable  at  non-zero  temperature  in  ID.  We  constructed 
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a  model  capturing  the  essential  properties  of  dipolar  fermions  in  optical  lattices.  We 
found  that  strong  anisotropic  interactions  suppress  unwanted  phase  disorder 
(fluctuations)  in  2D  to  pave  the  way  for  observing  these  interesting  excitations 
(Publication  13). 

•  Flat  bands  possess  macroscopic  level  degeneracy  and  play  a  crucial  role  in  many 
important  physical  phenomena.  However,  most  previous  lattice  models  for 
generating  flat  bands  involve  either  high  orbital  bands  or  high  order  tunnelings, 
which  are  generally  very  challenging  in  experiments.  We  showed  that  the 
combination  of  spin-orbit  coupling,  Zeeman  field  and  optical  lattice  potential  can 
yield  isolated  flat  ground  state  bands  where  topological  properties  may  originate 
from  the  spin-orbit  coupling  (Publication  14). 

•  Intriguing  many-body  states  are  possible  with  spin-orbit  coupling  but  implemented 
schemes  for  generating  the  necessary  fields  lead  to  unwanted  heating  and  losses 
near  Feshbach  resonances,  precisely  the  many-body  regime  of  interest.  We  found 
that  heating  and  losses  can  be  avoided  in  flat  bands  while  still  revealing  interesting 
fractionalized  many-body  states  (Publication  7). 

•  We  studied,  for  the  first  time,  the  BCS-BEC  crossover  in  a  spin-orbit  coupled 
degenerate  Fermi  gas.  We  showed  that  topological  excitations,  such  as  the  long- 
sought  Weyl  and  Majorana  fermions,  can  exist  in  such  systems.  We  found  various 
disorders,  including  the  low  dimensional  phase  disorders  (fluctuations),  do  not 
destroy  such  topological  states  (Publication  17,  20,  22). 

•  We  demonstrated  the  existence  of  a  new  route  for  realizing  the  long  sought  FFLO 
superfluids  using  spin-orbit  coupled  degenerate  Fermi  gases.  Such  FFLO  superfluids 
can  support  topological  excitations  such  as  Majorana  fermions  in  certain  parameter 
regimes.  We  further  showed  such  FFLO  phases  are  much  more  stable  than  the 
traditional  population  imbalanced  induced  FFLO  phases  in  low  dimension  (already 
studied  extensively  in  experiments),  where  the  finite  temperature  induced  phase 
disorder  (fluctuations)  are  crucial  (Publication  5,  9, 12). 

•  In  Collaboration  with  Peter  Engels  experimental  group,  we  observed,  for  the  first 
time,  the  long-sought  Zitterbewegung  phenomenon,  the  simultaneous  velocity  (thus 
position)  and  spin  oscillations,  for  neutral  atoms  using  spin-orbit-coupled  Bose- 
Einstein  condensate  (BEC)  through  sudden  quantum  quenches  of  the  system 
Hamiltonian.  Together  we  also  observed  a  moving  antiferromagnetic  order  using  a 
phase-winding  BEC.  Antiferromagnetic  order  has  been  extensively  studied  in  optical 
lattices  (one  major  goal  for  the  DARPA  optical  lattice  emulator  program),  and  only 
recently  some  signatures  were  observed  in  optical  lattice  experiments  (Hulet  group) 
(Publication  10, 11) 
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Emergent  Kinetics  and  Fractionally  Charged  Excitations  in  ID  Spin-Orbit  Coupled  Flat  Band 

Optical  Lattices, 

F.  Lin.  C.  Zhang,  and  V.  W.  Scarola 

Phys.  Rev.  Lett.  112, 110404  (2013)  arXiv:1308.5601 

Recent  ultracold  atomic  gas  experiments  implementing  synthetic  spin-orbit  coupling  allow 
access  to  flat  bands  that  emphasize  interactions.  In  our  work  we  modeled  spin-orbit  coupled 
fermions  in  a  one-dimensional  flat  band  optical  lattice.  We  introduced  an  effective  Luttinger 
liquid  theory  to  show  that  interactions  generate  collective  excitations  with  emergent 
kinetics  and  fractional  charge,  analogous  to  emergence  of  Fermi-liquid  properties  in  the 
fractional  quantum  Hall  regime.  Observation  of  these  excitations  would  provide  an 
important  platform  for  exploring  exotic  quantum  states  derived  solely  from  interactions. 
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Figure  Caption:  (Left)  Symbols  denote  the  man-body  energy  dispersion  computed  for  the 
Hubbard  interaction  in  a  flat  spin-orbit  band.  The  schematics  denote  a  real  space  depiction 
of  a  Wigner  crystal  ground  state  with  dispersive  collective  modes.  (Right)  Emergent  single¬ 
particle  energy  dispersion  versus  wavevector.  The  inset  depicts  the  emergent  hopping  that 
derives  purely  from  interactions. 
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Neel  temperature  and  thermodynamics  of  the  half-filled  3D  Hubbard  model  by 
Diagrammatic  Determinant  Monte  Carlo, 

E.  Kozik,  E.  Burovski,  V.W.  Scarola,  M.  Troyer 

Phys.  Rev.  B  87,  205102  (2013).  arXiv:  1212.3027 

Ultracold  atomic  gases  experiments  seeking  to  emulate  the  Fermi  Hubbard  model  are 
currently  attempting  to  reach  temperatures  low  enough  to  realize  spin  ordering.  In  this 
work  we  studied  the  thermodynamics  of  the  3D  Hubbard  model  at  half  filling  on  approach 
to  the  Neel  transition  by  means  of  large-scale  unbiased  Diagrammatic  Determinant  Monte 
Carlo  simulations.  We  obtain  the  transition  temperature  in  the  strongly  correlated  regime, 
as  well  as  temperature  dependence  of  energy,  entropy,  double  occupancy,  and  the  nearest- 
neighbor  spin  correlation  function.  Our  results  improve  the  accuracy  of  previous  unbiased 
studies  and  present  accurate  benchmarks  in  the  ongoing  effort  to  realize  the 
antiferromagnetic  state  of  matter  with  ultracold  atoms  in  optical  lattices. 
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Figure  Caption:  (Left)  Schematic  phase  diagram  showing  several  phases  of  high 
temperature  superconductors.  Arrows  indicating  current  experiments  and  our  theory  work 
show  the  operating  temperatures  of  optical  lattice  experiments  emulating  the  Hubbard 
model  and  our  calculations,  respectively.  (Right)  Nearest  neighbor  spin  correlation 
functions  and  the  double  occupancy  plotted  as  a  function  of  temperature  as  the  Neel  state  is 
approached  on  the  left  and  right  axes,  respectively. 
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Double  occupancy 


Observation  of  Zitterbewegung  in  a  spin-orbit-coupled  Bose-Einstein  condensate 

C.  Ou.  C.  Hamner,  M.  Gong.  C.  Zhang,  P.  Engels 

Phys.  Rev.  A,  88,  021604(R)  (2013).  (49  citations) 

Spin-orbit-coupled  ultracold  atoms  provide  an  intriguing  new  avenue  for  the  study  of  rich 
spin  dynamics  in  superfluids.  In  this  work,  we  observed  Zitterbewegung,  the  simultaneous 
velocity  (thus  position)  and  spin  oscillations,  of  neutral  atoms  between  two  spin-orbit- 
coupled  bands  in  a  Bose-Einstein  condensate  (BEC)  through  sudden  quantum  quenches  of 
the  Hamiltonian.  The  observed  Zitterbewegung  oscillations  are  perfect  on  a  short  time  scale 
but  gradually  damp  out  on  a  long  time  scale,  followed  by  sudden  and  strong  heating  of  the 
BEC.  As  an  application,  we  also  demonstrated  how  Zitterbewegung  oscillations  can  be 
exploited  to  populate  the  upper  spin-orbit  band  and  observe  a  subsequent  dipole  motion. 
Our  experimental  results  were  corroborated  by  a  theoretical  and  numerical  analysis  and 
showcased  the  great  flexibility  that  ultracold  atoms  provide  for  investigating  rich  spin 
dynamics  in  superfluids. 


Figure  Caption:  (Left)  Spin-orbit  coupled  band  structure.  The  sudden  change  of  the  detuning 
of  two  spins  induces  the  coupling  between  two  bands,  leading  to  Zitterbewegung 
oscillations.  (Middle)  The  experimental  observed  short  time  ZB  oscillation.  The  velocity  of 
atoms  oscillates  with  time.  (Right)  Decay  of  the  ZB  oscillation  in  a  long  time  scale  due  to 
many-body  interaction  between  atoms. 
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Route  to  Observable  Fulde-Ferrell-Larkin-Ovchinnikov  Phases  in  3D  Spin-Orbit  Coupled 

Degenerate  Fermi  Gases, 

Z.  Zheng,  M.  Gong.  X.  Zhou,  C.  Zhang,  and  G.-C.  Guo, 

Phys.  Rev.  A  87,  031602(R)  (2013).  arXiv:1208:2029  (37  citations) 


The  Fulde-Ferrell-Larkin-Ovchinnikov  (FFLO)  phase  was  first  predicted  in  two-dimensional 
superconductors  about  50  years  ago,  but  so  far  unambiguous  experimental  evidence  is  still 
lacking.  The  recently  experimentally  realized  spin-imbalanced  Fermi  gases  may  potentially 
unveil  this  elusive  state,  but  they  require  very  stringent  experimental  conditions.  In  this 
work,  we  showed  that  FFLO  phases  may  be  observed  even  in  a  three-dimensional  (3D) 
degenerate  Fermi  gas  with  spin-orbit  coupling  and  an  in-plane  Zeeman  field.  The  FFLO 
phase  is  driven  by  the  interplay  between  the  asymmetry  of  the  Fermi  surface  and  the 
superfluid  order,  instead  of  the  interplay  between  magnetic  and  superconducting  order  in 
solid  materials.  The  predicted  FFLO  phase  exists  in  a  giant  parameter  region,  possesses  a 
stable  long-range  superfluid  order  due  to  the  3D  geometry,  and  can  be  observed  with  an 
experimentally  already  achieved  temperature  (T  ~  0.05EF),  thus  opening  a  fascinating 
avenue  for  exploring  FFLO  physics. 
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Figure  Caption:  BCS-BEC  crossover  phase  diagrams  in  the  presence  of  SO  coupling  and  an 
in-plane  Zeeman  field,  (a)  Without  SO  coupling.  The  circle  symbol  represents  the  data  from 
the  quantum  Monte  Carlo  calculation,  (b)  and  (c)  with  SO  coupling  strength  ocKf  =  0.5Ef  and 
aKF  =  I.OEf.  (d)  In  the  unitary  regime. 
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Bose-Einstein  Condensates  in  Spin-Orbit  Coupled  Optical  Lattices:  Flat  Bands  and 

Superfluidity, 

Yongping  Zhang,  and  Chuanwei  Zhang 
Phys.  Rev.  A  87,  023611  (2013).  arxiv:1203.2389 


Spin-orbit  (SO)  coupled  superfluids  in  free  space  or  harmonic  traps  have  been 
extensively  studied.  However,  the  rich  physics  of  SO  coupled  BEC  in  optical  lattices  has  not 
been  explored.  In  this  project,  we  propose  an  experimentally  feasible  route  for  generating 
isolated  flat  bands  using  cold  atoms  in  SO  coupled  weak  optical  lattices.  Flat  bands  possess 
macroscopic  level  degeneracy  because  of  their  flat  energy  dispersion,  and  play  a  crucial  role 
in  many  important  physical  phenomena.  However,  most  previous  lattice  models  for 
generating  flat  bands  involve  either  high  orbital  bands  or  high  order  tunnelings,  which  are 
generally  very  challenging  in  experiments.  We  show  that  the  combination  of  SO  coupling, 
Zeeman  field  and  optical  lattice  potential  can  yield  isolated  flat  bands  where  topological 
properties  may  originate  from  the  SO  coupling. 


Our  proposed  SO  coupling  mechanism,  when  generalized  to  2D,  may  provide  an 
experimentally  feasible  route  for  generating  chiral  flat  bands  and  studying  relevant 
fractional  quantum  Hall  insulator  physics.  The  stable  superfluidity  in  the  whole  ground  state 
band  may  lead  to  other  interesting  phenomena  that  have  not  been  explored  in  regular 
optical  lattices,  such  as  dissipationless  Bloch  oscillation  of  BEC. 


Figure  Caption:  (Left)  The  experimental  scheme  for  generating  flat  band  in  spin-orbit  coupled 
optical  lattices.  (Right)  The  flat  band  energy  dispersion  (inset)  and  the  band  flatness  ratio  R  with 
increasing  Raman  coupling. 
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Percolation  Enhanced  Supersolids  in  the  Extended  Bose-Hubbard  Model, 
B.M.  Kemburi  and  V.  W.  Scarola 


Phys.  Rev.  B.  85,  020501  (2012)  arXiv:1110.2957 

In  this  work  we  theoretically  studied  the  stability  of  lattice  supersolid  states  in  the  extended 
Bose-Hubbard  model  with  bounded  spatial  disorder.  We  constructed  a  disorder  mean  field 
theory  and  compared  with  quantum  Monte  Carlo  calculations.  The  supersolid  survives  weak 
disorder  on  the  simple  cubic  lattice.  We  also  found  that  increasing  disorder  strength  can 
transform  a  lattice  solid  into  a  supersolid  as  it  tends  to  percolate  through  the  disorder 
landscape. 


Figure  Caption:  (Left)  Quantum  Monte  Carlo  results  for  a  single  disorder  configuration  for 
the  extended  Bose-Hubbard  model.  We  find  solid  order  among  all  sites  but  we  only  plot  a 
site  if  it  has  large  density  fluctuations.  The  supersolid  is  allowed  to  percolate  from  one  edge 
to  another  only  if  local  density  fluctuations  permit  local  superfluid  order  along  enough 
connected  bonds  to  cross  the  sample.  (Right)  The  top  and  bottom  panels  plot  order 
parameters  to  identify  the  superfluid  using  mean  field  theory  and  Quantum  Monte  Carlo. 

The  left  and  right  axis  plot  static  structure  factor  and  the  superfluid  density,  respectively,  for 
the  disordered  extended  Bose-Hubbard  model  on  a  cubic  lattice  as  a  function  of  the  disorder 
strength.  The  supersolid  arises  when  both  order  parameters  are  non-zero. 
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Enhancing  the  thermal  stability  of  entanglement  between  Majorana  fermions  with  dipoles  in 

optical  lattices, 

Fei  Lin  and  V.W.  Scarola 

Phys.  Rev.  Lett.  Ill,  220401  (2013)  arXiv:1210.0799 

Pairing  between  spinless  fermions  can  generate  Majorana  fermion  excitations.  Such 
excitations  may  exhibit  intriguing  properties  arising  from  non-local  entanglement,  including 
anyonic  braid  statistics,  teleportation,  and  enough  stability  to  encode  quantum  information. 
But  simple  models  indicate  that  non-local  entanglement  between  Majorana  fermions 
becomes  unstable  at  non-zero  temperatures.  We  addressed  this  issue  by  showing  that 
anisotropic  interactions  between  dipolar  fermions  in  optical  lattices  can  be  used  to 
significantly  enhance  thermal  stability.  We  constructed  a  model  of  oriented  dipolar 
fermions  in  a  square  optical  lattice.  We  found  that  domains  established  by  strong 
interactions  exhibit  enhanced  entanglement  between  Majorana  fermions  over  large 
distances  and  long  times  even  at  finite  temperatures. 


Figure  Caption:  (Left)  Schematic  of  oriented  fermionic  dipoles  in  an  optical  lattice.  Each  x- 
row  sets  up  a  paired  state  (a  Kitaev  chain)  with  Majorana  fermions  at  the  end  of  each  row. 
(Right)  The  fermion  parity  for  Kitaev  chains  plotted  as  a  function  of  a  perturbing  field  for 
two  different  temperatures  T.  The  high  temperature  plot  shows  the  stability  of  Majorana 
fermion  entanglement  against  thermal  fluctuations. 
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Searching  for  Majorana  Fermions  in  2D  Spin-orbit  Coupled  Fermi  Superfluids  at  Finite 

Temperature, 

Ming  Gong.  Gang  Chen.  Suotang  Jia,  Chuanwei  Zhang 

Phys.  Rev.  Lett.  109, 105302  (2012)  arXiv:1201.2238  (47  citations) 

Recent  experimental  breakthroughs  in  realizing  spin-orbit  (SO)  coupling  for  cold  atoms  has 
spurred  considerable  interest  in  the  physics  of  2D  SO  coupled  Fermi  superfluids,  especially 
topological  Majorana  fermions  (MFs)  which  were  predicted  to  exist  at  zero  temperature. 
However,  it  is  well  known  that  long-range  superfluid  order  is  destroyed  in  2D  by  the  phase 
fluctuation  at  finite  temperature  and  the  relevant  physics  is  the  Berezinskii-Kosterlitz- 
Thouless  transition.  In  this  project,  we  examined  finite  temperature  effects  on  SO  coupled 
Fermi  gases  and  showed  that  finite  temperature  is  indeed  necessary  for  the  observation  of 
MFs.  MFs  are  topologically  protected  by  a  quasiparticle  energy  gap  which  was  found  to  be 
much  larger  than  the  temperature.  The  restrictions  to  the  parameter  region  for  the 
observation  of  MFs  have  been  obtained. 
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BCS-BEC  Crossover  and  Topological  Phase  Transition  in  3D  Spin-Orbit  Coupled 

Degenerate  Fermi  Gases, 

M.  Gong.  S.  Tewari,  and  C.  Zhang, 

Phys.  Rev.  Lett.  107, 195303  (2011).  arXiv:1105.1796  (125  citations) 

Ultra-cold  degenerate  Fermi  gases  (DFG)  with  population  imbalance  (equivalent  to  an 
effective  Zeeman  field)  have  sparked  tremendous  recent  interest  in  both  theory  and 
experiment.  The  recent  experimental  realization  of  spin-orbit  coupling  (SOC)  opens  another 
new  avenue  for  the  study  of  DFG.  We  investigated,  for  the  first  time,  the  BCS-BEC  crossover 
physics  in  DFG  with  both  Zeeman  field  and  SOC.  The  crossover  physics  is  important  because 
the  s-wave  superfluid,  together  with  SOC  and  Zeeman  field,  may  yield  intriguing  chiral  p- 
wave  physics  with  non-trivial  statistical  properties.  For  example,  Majorana  fermions  and  the 
associated  non-Abelian  statistics,  the  critical  ingredients  for  performing  TQC,  may  be 
observable  only  in  the  crossover  region  where  the  superfluid  order  parameter  is  large  and 
thus  robust  against  finite  temperature  effects.  We  showed  that  the  superfluid  order 
parameter  destroyed  by  a  large  Zeeman  field  can  be  restored  by  SOC.  At  zero  temperature, 
we  found  that  there  is  a  topological  phase  transition  in  such  systems  from  a  non-topological 
superfluid  to  a  topological  superfluid  with  increasing  Zeeman  fields  or  by  tuning  the  s-wave 
scattering  length.  The  two  phases  are  separated  by  a  critical  point  which  itself  is  a  gapless 
superfluid.  We  showed  how  the  zero-temperature  topological  quantum  phase  transition  can 
be  probed  at  finite  temperatures  using  cold  atom  photoemission  spectroscopy  which  has 
already  been  realized  in  experiments. 


Figure  Caption:  (Left)  Change  of  the  order  parameter  in  BCS-BEC  crossover  with  spin-orbit 
coupling.  (Right)  The  Weyl  fermions  in  quasiparticle  excitations  and  the  corresponding 
topological  phase  transitions. 
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